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ABSTRACT 

We present an alternative explanation for the nature of turbulence in molecular clouds. 
Often associated with classical models of turbulence, we instead interpret the observed 
gas dynamics as random motions, induced when clumpy gas is subject to a shock. 
From simulations of shocks, we show that a supersonic velocity dispersion occurs in 
the shocked gas provided the initial distribution of gas is sufficiently non-uniform. We 
investigate the velocity size-scale relation a oc for simulations of clumpy and fractal 
gas, and show that clumpy shocks can produce realistic velocity size-scale relations 
with mean a ^ 0.5. For a fractal distribution, with a fractal dimension of 2.2 similar 
to what is observed in the ISM, we find a ex r^*^. The form of the velocity size-scale 
relation can be understood as due to mass loading, i.e. the post-shock velocity of the 
gas is determined by the amount of mass encountered as the gas enters the shock. We 
support this hypothesis with analytical calculations of the velocity dispersion relation 
for different initial distributions. A prediction of this model is that the line-of sight 
velocity dispersion should depend on the angle at which the shocked gas is viewed. 
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1 INTRODUCTION 

Molecular c louds are know n to ex h ibit supersonic chaotic 
dyna m ics ( Hever & Brunt 2004; Ossenkopf Mac Lowl 

■ | 2002l: iFalgarone Phillips 1990 : Perault et all Il985l: 

■ iHavashi et al.l Il989l : iLarsonl which are thought to 
I control star format ion and d etermine the p roperties of 

protostellar cores (|Mac Low Klesse^ 120041) . Although 
referred to as 'turbulence', the origin and nature of these 
motions are not fully understood. The most general def- 
inition of ISM turbulence is simply t hat the gas exhibits 
random motions on many scales (|Mac Low KlessenI 
120041 1. However there is a consistent correlation observed 
in molecular clo uds bet w een t he velocity dispersion and 
size scale ft he ILarsonl (Il98l|) relation), approx i mately 
a (X r °-^ (e.g. iMveral (Il98.1) : ISolomon et all (Il98/i) : iBruntI 



(|2003)). This has invoked many comparisons between 
interstellar turbulence and classical turbulence, e.g. Kol- 
morogov incompressible turbulence (|Passot et al.l Il988l : 
iFalgarone Philhpsl Il990h (a oc L^^)- Burger's shock 
dominated turbulence (jScalo et al.l Il998l) {a oc L^'^); 
and the She-Leveque model for incompressible turbu- 
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lence (IShe Leveaud Il994l : iBoldvrevI l2002^ (a oc L^'^^ 
JBoldvrev et al.ll2002M . 

The possible sources of turbulence can be summarised 
as follows: gravitational, magnetic or hydromagnetic insta- 
bilities; galactic rotation, through magneto-rotational in- 
stabilities, shocks in spiral arms or collisions of clouds on 
different epicyclic orbits; stellar feedback via supernovae, 
stellar winds and HII regions. Recent simulations have 
indicated that turbulence induced by a large scale driv- 
ing force (e.g. large scale flows from supernovae or galac- 
tic rotation) is more consistent with observ ed molecular 
cloud structures (|Bruntl l2003l : iKlessenlEool . Supernovae 
have been shown to produce sufficient energy to gener- 
ate t he velocity dispersions observed (Mac Low & KlesserJ 
120041) . However observations of turbulent velocities in re- 
gions which do not contain massive star formation suggests 
that ot her mechanisms, such as magnet o- r ot at ional inst a- 
bilities fPiontek k Ostriker 2005; Sellwood k Balbus 1999) 
and colliding flows ( Ballesteros-Pare des et a l. 1999) must 
also be important. Interestingly, recent observations have 
suggested that the elongations of molecular clouds are more 
compatib le with galactic ro tation models rather than stellar 
feedback (iKoda et al.ll200(tl . 

Galactic disk simulations have investigated gravity 
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driven turbulence dWac a et al.l 1200211. stellar fee dback 
(IWada Nori^ l200ll: Ide Avillez Breitschwerdtl Eooi 
IPib et ani200(t) and the influ ence of spiral density waves on 
ISM dynamics (|Dobbs et al.l Eo06V Analytical results also 
indicate that vorticity is generated in centrally condensed 
clouds subject to galactic shocks (Kornreich & Scalo 2000), 
and the induced velocities follow the observed velocity size- 
scale relation. Previous numerical work on colliding flows 
showed that density and velocity perturbation s occur even in 
uniform flows subject to cooling instabilities (H eitsch et al.l 
I2OO 5V although a velocity length scale correlation was not 
investigated. Simulations of clumpy flows have also indicated 
that a Salpeter type clump mass spectrum can be repro- 
duced JClark ^, Bonnell 2006). 

Spiral shocks have also be en proposed to explain the dy- 
nami c s of molecular clouds (^Bonnell et al.l200^ : lzhang et al.l 
bOOlh . lBonnell et al.l()200rf) model giant molecular cloud for- 
mation as gas passes through a clumpy spiral shock. The dy- 
namics of the molecular clouds are determined on all scales 
simultaneously as the clouds form and the induced velocity 
dispersion size scale relation is consistent with observations. 
This can account for the observed velocity dispersions that 
are found even in regions devoid of massive stars. Further- 
more, there is no need for a continuous driving mechanism 
as the time for the decay of these velocities is proposed to 
be of similar magnitude to the cloud lifetime. 

In this paper we investigate the velocity size relation in 
shock tests with uniform, clumpy and fractal distributions 
of isothermal gas. The clumpy and fractal distributions are 
chosen to r e flect the highly structure d nature of the ISM 
(ICo^ l2004 lElme green Scald l2004l : IPickev SsT, LockmanI 
ll99Gl : IPerault et al. 1985). We concentrate on modelling the 
passage of gas through a spiral arm by using a linear sinu- 
soidal potential, although these results would apply gener- 
ally to shocks between colliding flows. We show that random 
velocities induced in non-uniform shocks display a velocity 
size relation similar to those observed and provide simple 
analytical analysis alongside the results of our simulations. 
Thus the 'turbulence' in our results describes random mo- 
tions of the gas and does not correspond to any theories of 
classical turbulence. 



Test Distribution L F (%) M 



ST 


Uniform 


4 


100 


10 


ST 


Clumps (rcz =0.1) 


4 


27 


20 


ST 


Clumps (rcz=0.2) 


4 


50 


20 


ST 


Uniform 


4 


100 


20 


ST 


Clumps (rci^O.l) 


4 


27 


40 


ST 


Clumps (rcz=0.2) 


4 


50 


40 


SP 


Uniform 


3 


100 


30 


SP 


Clumps (rc^ =0.1) 


3 


12 


30 


SP 


Clumps (rcz=0.2) 


3 


23 


30 


SP 


Clumps (rc;=0.4, 0.1, 0.04) 


3 


32 


30 


SP 


2.2 D Fractal 


3 


10 


30 


SP 


2.7 D Fractal 


3 


23 


30 



Table 1. Table showing the different runs performed and the 
distributions of gas used for initial conditions. ST and SP refer 
to shock tube and sinusoidal potential tests, L is the maximum 
initial length scale, F the volume filling factor and M. the Mach 
number of the shock. All simulations used 2 x 10^ particles except 
for the fractal distributions which required either 2 x 10^ or 3 x 10^ 
particles. 



characterized by the Mach number of the shock and the ini- 
tial density distribution of the gas. In all calculations, the 
particles are allocated velocities in the x direction only and, 
except for the obhque shocks (Section 3.2), the gas shocks 
in the yz plane. 

The parameters used in the different simulations are 
shown in Table 1. We investigated 4 distributions of gas - 
uniform, homogeneous spherical clumps in pressure equilib- 
rium, spherical clumps of different radii/density and frac- 
tal distributions. The filling factors for each distribution 
are calculated by overlaying a 3D grid on each distribution 
and determining the porosity. The filling factor is given by 
F — Nfuii/Ntotai where Nfuii is the number of cells con- 
taining at least 1 particle and Ntotai the total number of 
cells. We take a 32^ grid, so the mesh resolution is equiva- 
lent to a couple of smoothing lengths for the uniform case. 
This ensures that each cell contains at least 1 particle for 
the uniform distribution so the filling factor is 100%. The 
maximum initial scale length corresponds to the range of x 
values for the particles in the initial distribution. 



2 CALCULATIONS 

We use the 3D smoothed partic le hydrody namics (SPH) 
code based on the version by Benz (lBenzll990h . The smooth- 
ing length is allowed to vary with space and time, with the 
constraint that the typical number of neighbours for each 
particle is kept near N neigh ~ 50. Artificial viscosity is in- 
cluded with either the standard parameters q = 1 and [3 — 2 
fMonag han Lattanziolll985l : lMonaghanlll992^ or a = 2 
and (3 = A. 



2.1 Initial conditions 

We investigate uniform, clumpy and multi-scale shocks, 
starting with 3D shock tube tests first before considering 
gas subject to an external potential. In all calculations, the 
gas is isothermal, non self-gravitating and there are no mag- 
netic fields. These calculations are dimensionless, and are 



2.1.1 Shock tube test 

We first perform 3D shock tube tests with initial distribu- 
tions of uniform and clumpy gas to model colliding flows. 
For both distributions, particles are placed within a cuboid 
of dimensions — 2<x<2, —l<y<l and — 1 < 2; < 1. To 
produce a clumpy shock we distribute the particles in uni- 
form density spheres within these length scales. The clumps 
are a constant temperature, conflned by either an external 
pressure fleld or a hotter diffuse phase. The clumps are ini- 
tially allowed to settle into equilibrium before the simula- 
tion is carried out. To produce clumps of different radii, the 
external pressure (or pressure of the diffuse phase) can be 
increased or decreased. We then assign particles a velocity of 
vq = 10 Cs for X ^ and = — 10 Cs for x > 0. For the uni- 
form distribution, this produces two approximately Mach 
10 shocks, and multiple shocks of up to Mach 20 for the 
clumpy distributions. These calculations were also repeated 
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using twice the initial velocities and all tests use 2 x 10^ 
particles. 

2.1.2 Sinusoidal potential 

We mainly consider shocks in gas subject to an external 
potential. In this case the gas self shocks, similar to gas ex- 
periencing a stellar potential in a spiral galaxy. A velocity 
dispersion relation similar to the observed a oc r°"^ law has 
been sho wn to develop as gas passes through a clumpy spi- 
ral shock (|Bonnell et al.l2QQ(t ) . Here we examine a simplified 
setup, where we can investigate the effect of the shock dy- 
namics on the initial gas distribution. Instead of the spiral 
potential we use a ID sinusoidal potential of the form 

= Acos{k{x + B)). (1) 

where k is the wavenumber and B a length parameter to 
determine the location of the minimum. This is equivalent to 
the linear passage of gas through sinusoidal spiral arms. The 
dimensionless velocity acquired by gas falling from the peak 
to the base of the potential is Vpot = V2 x A. We tried many 
different potentials, varying A and k, and only applying the 
potential once the gas has passed a minimum. However we 
found that the results presented here are largely independent 
of the exact nature of the potential. The structure of the 
shock is similar for different potentials for a given initial 
distribution. The relative strength of the shock determines 
the magnitude of the velocity dispersion, whilst the initial 
distribution determines the velocity size scaling law. For the 
simulations presented here, we took /c = 7r/4, A = 100 and 
B = 2 to produce a minimum at 2 and maxmima at -2 and 
6. 

We allocate particles a velocity of 50 Cs in the x di- 
rection and zero velocity in the y and z directions, which 
for the simulations described here, leads to a shock of Mach 
number ^ 30. We set up a distribution of spherical clumps 
in pressure equilibrium in the same way as described for the 
shock tube tests. Where a hot diffuse phase is used to supply 
an external pressure, the hot phase is distributed with the 
same number of particles (2x10^), but 1/10 of the mass of 
the cold phase. 

We also test a distribution with clumps of different size- 
scales and densities, giving structure on a range of scales. 
The clumps have initial diameters of 0.4, 0.1 and 0.04, with 
clumps of smaller diameter placed inside larger clumps. For 
the uniform and clumpy distributions, particles are posi- 
tioned within a cuboid of dimensions —1.5 < x < 1.5, 
— 1 < y < 1 and —l<z<l. We also investigate shocks 
with an initial fractal distribution, following the method de- 
scribed in Elmeereen (1997) to generate fractals. The algo- 
rithm includes 3 parameters, an intrinsic length scale L, the 
number of hierarchical levels, H, and the number of points in 
each level, N. The dimension of the fractal is D = logN/logL 
and the number of points is . We generate a 2.2 D fractal, 
with L — 2.1, N — 5 and H — 9 requiring ^ 2 million points, 
and a 2.7 D fractal, with L — 2.5, = 12 and H — Q requir- 
ing ^ 3 million points. We then scale the x,y, z coordinates 
(equally) of each fractal to fit inside a cube of dimensions 
-1.5 <x < 1.5, -1.5 < y < 1.5 and -1.5 < z < 1.5. Obser- 
vations estimate the inte rstellar fractal dimension as D=2.3 
(jElmegreen Falgarond [l996'). 

Since for the distributions with fractals or different 



size clumps the gas exhibits different densities and pres- 
sures, constant pressure boundaries, or an intervening dif- 
fuse phase, are no longer appropriate. Instead we apply a 
pressure switch in order that only gas in the shocked region 
is subject to pressure forces. The gas experiences pressure 
only when div{v) ^ i.e. compression of the gas is occurring. 
This enables structure on all scales to be maintained in the 
gas distribution before gas reaches the shock. Tests for the 
uniform density clumps showed this method produced sim- 
ilar results compared to when constant pressure boundaries 
were applied. 



3 RESULTS 

Column density plots for different simulations are shown in 
Fig. 1,2,3 (all for the sinusoidal potential tests). The uniform 
shock in Fig. 1 shows a smooth shocked region of approxi- 
mately constant density and width. By contrast the clumpy 
shock (Fig. 1, middle) shows a much broader shocked re- 
gion of non- uniform density. The shock contains more struc- 
ture and appears more similar to simulations of turbulence. 
In Fig. 1 (middle), an external pressure field is applied to 
maintain the clumps in pressure equilibrium. Fig. 1 (lower) 
shows a shock for similar size clumps, where the clumps are 
instead surrounded by hotter gas. In this case, the hot gas 
takes the same numbers of particles as the cold gas, but 1/10 
of the mass. The ratio of the densities of cold to hot gas is 
^ 30, which is similar to the ratio of densities of the cold 
and warm neutral components of the ISM (fCoxii200^ . The 
structure of the cold gas in the shock is very similar whether 
the clumps are in equilibrium from external pressure bound- 
aries, or hot gas, indicating that the hot gas has little effect 
on the gas dynamics. 

The distribution of different size clumps shows similar 
morphology, although more smaller scale structure is appar- 
ent in the shocked gas (Fig. 2). The shocked gas of the frac- 
tal distribution (Fig. 3) shows more filamentary structure 
compared to the clumpy distributions. 

3.1 Velocity dispersion 

We now calculate the ID velocity dispersion of the post- 
shock gas. We only consider the Vx velocities, which corre- 
sponds to the direction of motion of the initial gas, since 
the velocity dispersion in the y and z directions are always 
subsonic. For a given size, we average the velocity disper- 
sion over numerous regions of that size scale. The regions 
are 3 D and chosen to centre on the densest particles in the 
shock. Only particles with densities greater than the max- 
imum pre-shock density are considered for calculating the 
velocity dispersion, thus ensuring we only include gas in the 
shock. (We find that even for a uniform shock, including the 
pre and/or post shock gas will produce a Larson type ve- 
locity dispersion size-scale relation, theoretically and from 
numerical results.) We repeat this process for regions of dif- 
ferent size-scale to determine the dependence of the induced 
velocity dispersion on the size-scale. 

For the shock tube tests, we initially found that the ve- 
locity dispersion was supersonic, even for the uniform shock. 
This is due to the inherently clumpy nature of SPH which in- 
troduces error when calculating the velocity dispersion. We 
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Figure 1. The column density distribution is shown for hnear 
shock tests where gas is passed through a sinusoidal potential. The 
initial distribution is uniform (top) and clumpy (middle, bottom, 
where the clump radius is 0.1). The middle panels show the case 
where the clumps are confined by constant pressure boundaries, 
whereas for the bottom panels, the clumps are in equilibrium with 
a diffuse phase of gas. The left panels (2.8 < x < 4.8, — 1 < y < 1) 
show the gas as it begins to shock (t=0.15) while the right panels 
(3.2 < X < 5.2, — 1 < 2/ < 1) show the gas when the shock is fully 
developed (t=0.25). The minimum and maximum of the potential 
lie dX X = 2 and x = Q respectively. The column density ranges 
form 90 to 9 X 10^ in units of particles per unit area and the 
colour scale is logarithmic. 

therefore increased the viscosity parameters to a — 2 and 
[3 — A. This lowered the values of the velocity dispersion for 
both the uniform Mach 10 and Mach 20 shocks, although 
the velocity dispersion is still supersonic for the Mach 20 
shock. The results presented for the shock tube tests use the 
higher viscosity parameters, whilst for the sinusoidal poten- 
tial tests, there is less noise and the standard parameters 
a — 1 and (3 — 2 are used. In all cases the higher viscos- 
ity parameters had little effect on the velocity dispersions 
for the clumpy shocks. Alternatively, the velocity dispersion 
can be determined from the SPH smoothed velocities: 

^'=^' + ET^^yWy. (2) 

Since these velocities are smoothed over the neighbouring 
particles, the velocity dispersion produces less noise. Conse- 




Figure 2. The column density distribution is shown for a 
shock test where the initial distribution consists of different radii 
clumps. The gas shocks as it passes through a sinusoidal potential. 
The left panel ( — 1 < x < 1, —1 < y < 1) shows the initial gas dis- 
tribution (t=0) while the right panel (3.2 < x < 5.2, -1 < y < 1) 
shows a stage (t=0.25) during the shock. The scaling is the same 
as Fig. 1, and the minimum and maximum of the potential lie at 
X = 2 and x = 6. 




Figure 3. The column density distribution is shown for a shock 
test where gas distributed according to a 2. 2D fractal is passed 
through a sinusoidal potential. The left panel ( — 1.5 < x < 1.5, 
— 1.5 < y < 1.5) shows the initial gas distribution (t=0) while the 
right panel (3.2 < x < 5.2, —l<y< 1) shows a stage (t=0.25) 
during the shock. The scaling is the same as Fig. 1. 



quently the velocity dispersion is lower (and in all our results 
subsonic) for the uniform shocks. Except for Fig. 6 though, 
we use the SPH velocities, since these are the velocities pro- 
duced by the code. 

We plot the velocity dispersion size-scale relation for 
each of the simulations in Fig. 4, 5, 6, 7 and 8. Also shown 
is the a oc r°'^ relation, coinciding with most observational 
results. In Fig. 4 and 5, we show the velocity size-scale re- 
lation for uniform and clumpy initial gas distributions, for 
the shock tube test and the sinusoidal potential. For the si- 
nusoidal potential, the velocity dispersion for the uniform 
shock dispersion remains flat and subsonic (the sound speed 
in all simulations is 0.3). This is as expected, since for a uni- 
form shock, the velocity of the shocked gas should have zero 
velocity dispersion. Again for the shock tube tests, the ve- 
locity dispersions for the uniform shocks are relatively flat. 
By contrast the clumpy shocks show an increasing veloc- 
ity dispersion with size-scale. The velocity of gas in the 
shock depends on the amount of mass it has encountered 
(Section 3.3). For the clumpy shock, gas entering the shock 
will encounter different amounts of mass (e.g. where gas ap- 
proaches another clump, or alternatively a relatively empty 
area) and a range of velocities are exhibited by the shocked 
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Figure 4. The one-dimensional velocity dispersion of the post-shock gas is plotted against size-scale for the 3D shock tube test. The 
figures show a A4=20 shock (left) and A4=40 shock (right). The initial distributions prior to the shock are uniform (dotted line), clumps 
of radii 0.1 (solid line) and clumps of radii 0.2 (dash-dot line). External pressure boundaries are applied to keep the clumps in equilibrium. 
The dashed line represents the observed a oc r^'^ relation and the sound speed is 0.3. The initial distribution has a maximum spatial 
extent of 4, and the viscosity parameters a = 2 and [3 = 4 are used. 




Figure 5. These panels show the one-dimensional velocity dispersion of the post-shock gas plotted against size-scale for the sinusoidal 
potential. The left panel uses the same initial distributions (with the same key) prior to the shock as Fig. 4, although the maximum 
initial length scale is 3 for these distributions. The right panel compares the velocity dispersion when external pressure boundaries are 
applied (solid line) and a diffuse intervening medium is used (dot dash line). For the right panel, the clump radius is 0.1, so these results 
correspond to the column density images in Fig. 1 (middle and bottom). The dashed line shows a oc r^-^ and the sound speed is 0.3. 




0.1 , /, ,<r; J 

0.01 0.10 1.00 

Size 

Figure 6. The ID velocity dispersion relation for is plotted for post- 
shock gas using the actual and smoothed velocities. The results are 
shown for the r=0.1 clumpy distribution (solid line, actual veloci- 
ties, dot-dash line, smoothed velocities) and the uniform distribu- 
tion (dotted line, actual velocities, dot-dot-dot dash line, smoothed 
velocities), where the gas passes through the sinusoidal potential. 
The dashed line shows a oc r^-^ and the sound speed is 0.3. 




0.01 0.10 1.00 

Size 

Figure 7. The ID velocity dispersion dependence on size-scale is 
plotted for post-shock gas using an initial distribution of different 
size clumps. The maximum initial length scale is 3 and the gas shocks 
when passed through a sinusoidal potential. The dashed line show 
cr oc r^-^ and the sound speed is 0.3. 
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0.1 J 

0.01 0.10 1.00 

Size 

Figure 8. The ID velocity dispersion dependence on size-scale 
is plotted for post-shock gas where the initial distributions are a 
2. 2D (solid) and 2.7 D (dotted) fractal. The gas is passed through 
a sinusoidal potential and the maximum spatial extent of the 
initial fractals is 3. The dashed line show a oc r^-^ and the sound 
speed is 0.3. 

gas. This range of velocities increases as the size-scale of 
the region increases. At some size-scale, the region of gas 
will contain the full range of structure inherent in the initial 
distribution. The velocity size-scale relation then remains 
relatively flat for any further increase in size-scale. 

Fig. 5 also compares the velocity dispersion relation 
when constant pressure boundaries are applied compared 
to using a diffuse hot medium. The average gradient from 
the 2 slopes is very similar. Although there is some differ- 
ence in the form of the velocity dispersion relation, the dy- 
namics appear to be dominated by the cold gas. Fig. 6 com- 
pares the velocity dispersion calculated using the actual and 
smoothed velocities, for a uniform and clumpy shock. Using 
the smoothed velocities reduces the velocity dispersion in 
the uniform shock by a factor of around 2. The velocity 
dispersion is reduced at small scales for the clumpy shock 
leading to a slight increase in the gradient. 

In reality, the ISM has structure on many different 
length scales. To explore how this affects the resulting ve- 
locity dispersion, we have run simulations with a range of 
clump sizes and with initially fractal distributions. Fig. 7 
and 8 also show an increasing velocity size-scale relation for 
initial gas distributions exhibiting structure on a range of 
scales. Fig. 7 shows the velocity size-scale dependence where 
the initial distribution consists of different sized clumps, and 
Fig. 8 the initial fractal distributions. The velocity disper- 
sion extends to smaller scales due to the presence of struc- 
ture initially over these scales. Again, using the smoothed 
velocities slightly increases the gradient of the velocity dis- 
persion against size-scale. The velocity size-scale relation is 
similar to the observed relation a oc r^'^ for most of the re- 
sults in Fig. 4, 5, 6, 7 and 8 corresponding to non-uniform 
initial distributions. The exponent a, where a oc r"", varies 
from approximately 0.29 for the 2.7D fractal distribution to 
0.75 for the distribution with different clump sizes. With the 
exception of the a — 0.75 result, our result s lie within the 
observed range of values e.g. a — 0.5 =b 0.05 ([Solomon et al.l 
11987*; 'Dame et al.""l986^, a = 0.4 d i 0.1, a = .65 ± 0.08 
(Fuller & Mvers 1992), a 0.2, 0.5 (iGoodman et al. 1998), 
representing a range of size scales from cloud cores to molec- 
ular clouds. 



Generally the distributions with higher filling factors 
produce shallower gradients in the velocity dispersion rela- 
tion. As the filling factor increases, the distributions and 
subsequent velocity dispersions tend towards those of uni- 
form gas, and subsequently the maximum velocity dispersion 
is less supersonic. The exception appears to be the distribu- 
tion of different sized clumps, where the velocity size scale 
relation is somewhat steeper, probably because the filling 
factor of the smaller clumps is lower, and the smaller clumps 
contain most of the mass. 

3.2 Oblique shocks 

Due to the geometry of the shock, the velocity dispersion 
relations shown in the previous section have only used the 
Vx component of the velocity. We investigate a more general 
case by modifying Equation (1) to include a dependence on 
y. This produces an oblique shock, more similar to spiral 
shocks. The potential is still sinusoidal, but the minima of 
the potential now lie in planes inclined by an angle to the 
yz. This is equivalent to an inclination of 90^^ — 6 with the 
initial flow of the gas. We perform a shock test with this 
modified potential using the clumpy (rcz =0.1) initial distri- 
bution. The gas column density for the shock is shown in 
Fig. 9. The shock is found to induce a velocity dispersion in 
both Vx and Vy although the Vz dispersion is still subsonic. 
In Fig. 10, the velocity dispersion is displayed for Vx and Vy, 
where the shock is inclined at 45° and 30°. The magnitude 
of the Vx and Vy dispersions is then proportional to the com- 
ponent of the shock front perpendicular to these directions, 
i.e. 

ax _ max {\a\ cos 0,Cs) . . 

ay max {\a\ sin 0,Cs) 

for a given size scale. The term Cs is an estimate of the min- 
imum value of the dispersion, along the line of sight parallel 
to the shock (though as shown next in Fig. 11, the minimum 
velocity dispersion in these simulations is slightly lower than 

Cs). 

In Fig. 11, we display the velocity dispersion along the 
line of sight for different viewing angles, again for the 30° 
and 45° oblique shocks. The maximum velocity dispersion is 
chosen at each viewing angle, corresponding to a size scale 
of ~ 1 (c.f. Fig. 10). The angle is measured anticlockwise 
from the x > axis, with 10° increments. The peaks then 
correspond to the line of sight perpendicular to the shock, 
whilst the minima occur when the line of sight is parallel to 
the shock. The difference in magnitude of the peaks occurs 
since the amplitude of the 45° shock is stronger. From these 
results, we would expect some anisotropy in the magnitude 
of the velocity dispersion of gas in molecular clouds, corre- 
sponding to the geometry of shocks in the gas. For classical 
models of turbulence, the magnitude of the velocity disper- 
sion is unlikely to show any preferred direction. 

3.3 Mass loading 

We now examine physically why an increasing velocity size- 
scale relation emerges in our models. For the clumpy shocks, 
the gas exhibits a range of densities and velocities across the 
region of shocked gas. The post-shock velocity of a small 
parcel of gas depends on the amount of mass it encounters 
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Figure 9. The column density distribution is shown for an 
obUque shock where the shock front is incUned at 45^ to the y 
axis. The gas initially has a clumpy (rci=0.1) distribution and is 
passed through a sinusoidal potential. The figure (2.6 < x < 4.6, 
—0.8 < y < 1.8) shows a stage (t=0.2) during the shock, and the 
shock in this image lies half way between a potential minimum 
and the following maximum. The scaling is the same as Figure 1. 
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Figure 11. The maximum velocity dispersion along the line sight 
is plotted against viewing angle cf) for the 45° (solid) and 30° 
(dotted) oblique shocks. The angle (j) is measured anticlockwise 
from the x > axis. 
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Figure 10. The Vx (solid) and Vy (dotted) velocity dispersions 
are plotted against size-scale for post-shock gas where the shock is 
inclined at 45° (top) and 30° (bottom) to the y axis. For = 45°, 
(Jx and ay are almost identical. The clumpy (rcz =0.1) distribu- 
tion is used, with a maximum initial length scale of 3. The gas 
shocks when passed through a sinusoidal potential. The dashed 
line shows a oc r^-^ and the sound speed is 0.3. 

during the shock ('mass loading'). If proportionate amounts 
of gas enter the shock, conservation of momentum deter- 
mines that the velocity of the gas in the rest frame of the 
shock will be small. However, if gas entering the shock en- 
counters only a small amount of material, it's velocity will 
be less affected and remain of higher magnitude. 

Within a region of size-scale less than the structures in 
the gas, gas in that region will encounter a similar column 
density in the shock. Therefore the gas will exhibit a low 



velocity dispersion. However as the size-scale of the region 
increases, the region will include different structures and gas 
of different densities. Therefore different parcels of gas will 
encounter different amounts of mass and exhibit different 
post-shock velocities. Thus over a larger region a higher ve- 
locity dispersion occurs in the gas. 

To test the hypothesis of mass loading, we use a semi- 
analytical approach to calculate the velocity dispersion for 
the fractal and clumpy distributions. We set up the initial 
distribution of particles and position a grid across the plane 
in which the shock will occur, at the centre of the distribu- 
tion. We then calculate the mass of particles located within 
a distance / either side of the grid for each grid cell. By 
applying the conservation of momentum, we determine the 
expected velocity of the gas in each grid cell, assuming that 
all the gas within the width / of the centre is compressed 
into the shock. We took /=0.2 for the clumpy shock and 
/=0.6 for the fractal distribution where more of the gas ends 
up in the shock (this parameter does not change the results 
providing sufficient material is included to accurately rep- 
resent the gas distribution). We then calculate the velocity 
dispersion as described in Section 3.1. 

In Fig. 13 we display the velocity size-scale relation for 
the clumpy distribution (clumps of radius 0.1) and the 2 
fractal distributions, from simulations and the correspond- 
ing semi-analytical tests. The results for the simulations 
are all taken from the sinusoidal potential tests. For the 
clumpy distribution the slope determined from mass load- 
ing is somewhat steeper than the results from the simulation. 
The shape and gradient for the 2.2 D fractal shows a strong 
correlation between the simulation and the analytical result. 
For the 2.7D fractal, the semi- analytical velocity dispersion 
size relation is shallower than that determined by the sim- 
ulation, although both show a shallower gradient compared 
to the other 2 distributions. The analytical method can be 
repeated by setting up distributions with different random 
seeds to show the degree of scatter in the expected gradi- 
ents for the distributions. The clumpy shock, which should 
give very similar distributions regardless of the initial seed, 
produced a consistent slope of a = 0.48 ± 0.03 compared 
to a ~ 0.43 for the simulation. The fractal distributions 
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Figure 12. The dependence of the one-dimensional velocity dis- 
persion of the post-shock gas on size-scale is plotted from our 
simulations (solid) (as in Fig. 5 and 8) and the corresponding 
semi-analytical result (dotted). The semi-analytical result is cal- 
culated assuming the post-shock velocities are dependent on mass 
loading. Three different distributions are shown: clumpy (clump 
radii 0.1) (top); 2. 2D fractal (middle) and 2.7D fractal (bottom). 
The dashed line shows a oc r^-^ . 



showed a much greater variation in slope, the 2. 2D fractal 
giving a — 0.43 ±0.16 compared to a ~ 0.4 for the simula- 
tion. For the 2.7D fractal, we found a — 0.25 it 0.15 for the 
analytical result, so the line shown on Fig. 12 is at the lowest 
extent of this range, whilst a ~ 0.3 for the simulation. 



4 ANALYTICAL MODELS 
4.1 Collision of two clumps 

The previous results represented plausible mass distribu- 
tions for the ISM, although consequently we were only able 



to calculate the velocity size-scale relation numerically. Here 
we consider the case of 2 clumps colliding and calculate 
directly the resulting velocity dispersion. This distribution 
represents the 'least uniform' distribution, so the most re- 
moved scenario from a uniform shock. The clumps are as- 
sumed to be uniform density spheres and are travelling with 
equal and opposite velocities. The clumps are offset from 
each other and so do not collide head on. The parameters 
describing the collision are shown in Fig. 13, where b is the 
impact parameter of the 2 clumps. In the following descrip- 
tion, we choose a Cartesian axis centred on the point of max- 
imum overlap between the 2 spheres, such that the centre of 
each sphere is situated at 2/ = ±6/2. In these coordinates, 
the collision of the clumps is symmetrical about z = 0. We 
take = 1, and the column density is scaled so that maxi- 
mum column density is Umax = 1 along the diameter of each 
sphere. We first assume that the radii of the clumps are the 
same, ri — 1 and r2 = 1. 

To calculate the velocity dispersion, we assume that all 
the gas in the 2 clumps is compressed when they collide. We 
calculate the post-shock velocity of the material by deter- 
mining, for each clump, the mass per unit area perpendic- 
ular to the shock (SI and E2 in Fig. 13). The calculation 
of the velocity dispersion is centred along the column which 
contains the most mass. In Fig. 13, this corresponds to the 
column of gas at ^ = 0, z = (i.e. the midpoint of h). 

In Fig. 14 we plot the column density for each clump and 
the post-shock velocity, for a cross-section of the collision 
along the y axis. By taking a cross-section, the problem is 
reduced to 1 dimension, and the mass for each sphere can 
be determined as a function of length scale. Using the idea 
of mass loading we assume that all gas is compressed into 
the xy plane. We calculate the expected velocity of gas in 
the shock by applying the conservation of momentum, i.e. 



(El - T.2)vQ 
(S1 + S2) ■ 



(4) 



In Fig. 14, b — 0.5 and the 2 spheres are overlapping by a 
distance of 1.5. At y = there is an equal contribution in 
mass from each sphere, and so the velocity in the shock will 
be 0. The velocity increases to ±1 for gas which does not 
encounter any material from the other clump. 

We extend this idea to calculate the velocity across a 
2D shock from the collision of the clumps. We again assume 
that the gas is compressed in the yz plane and calculate 
the mass density entering the shock from each sphere. We 
then determine the expected post-shock velocity of the gas, 
again over the yz plane. We calculate the velocity dispersion 
over a disk centred on the co-ordinates ^ = 0, z = 0. The 
velocity dispersion is determined over disks of increasing ra- 
dius, to produce a velocity size-scale relation. We plot the 
mass- weighted velocity dispersion in Fig. 15, with collisions 
of different impact parameter. 

For the simple case of 2 clumps colliding, the velocity 
dispersion increases with size scale approximately as a power 
law, independent of the impact parameter. Increasing the 
impact parameter transposes the power law to smaller size 
scales, since the size of the shocked region decreases. For 
size-scales > r — 6/2, the velocity dispersion includes gas 
which does not enter the shock and the velocity size-scale 
relation flattens. The power law for the shocked gas is some- 
what steeper than the a oc r°'^ relation observed, rather 
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Figure 13. Diagram showing a collision with impact parameter b 
between 2 spheres. The column densities El and E2 perpendicular 
to the shock are calculated to determine the expected velocity of 
gas in the shock. 




Figure 14. The velocity (solid line) is plotted for a cross-section 
of post-shock material when 2 spheres collide. The cross-section 
corresponds to taking a cut where z = 0. The dotted and dashed 
lines show the column density for each sphere along the cross- 
section. 
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Figure 15. The mass- weighted velocity dispersion predicted for 
the collision of 2 spherical clumps, each of radii 1 and initially 
travelling at vq = 1. The impact parameters of the collision are 
b = 0.1 (dotted), b = 0.5 (solid), b = 1 (dot-dash) and b = 1.75 
(long dash) whilst the short dash line represents a oc r^-^ . 
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Figure 16. The mass-weighted velocity dispersion predicted for 
the collision of 2 spherical clumps travelling with initial velocity 
vq = 1. The radius of the first clump is set to 1, whilst the radius 
of the second clump is 1 (solid), 0.75 (long dashed), 0.5 (dot-dash) 
and 0.25 (dotted). The impact parameter for all the collisions is 
6 = 1. The short dash line shows cr oc r*^"^ . 

a (X . However, the power law exponent may be expected 
to decrease when considering larger regions and structure on 
multiple scales composed of many clumps (Section 4.2). 

This power law is also independent of the comparative 
radii of the 2 spheres. Fig. 16 shows the effect of varying 
the radius of the second clump (r2). Similarly to varying 
the impact parameter, the velocity size-scale relation shifts 
to smaller scales as r2 decreases, and levels off at size-scales 
comparable to r2. Scaling both clump radii up or down will 
extend the velocity size-scale relation to larger or smaller 
size-scales. For example, if the radius decreases by a factor of 
10, the column density and therefore mass- weighted velocity 
dispersion will also decrease by the same factor. 

Assuming that the clumps have velocities vq = 1, the ve- 
locity dispersion for an impact parameter of 6 ^ 0.5 reaches 
a maximum of ~ 0.3. Since the initial velocity v can be scaled 
up or down, the maximum velocity dispersion for an initial 
velocity of vqCs will be approximately 0.3 vqCs (e.g. 3 Cs if 
M — 20) . When the impact parameter is small, the velocity 
dispersion is unlikely to be supersonic (e.g. the maximum 
is ~ 0.5 Cs when 6 = 0.1 and M = 20). This is expected, 
since \i h — the clumps collide head on and the velocity 
dispersion is everywhere. 

4.2 Multiple collisions of clumps 

The velocity size-scale relation obtained for 2 clumps collid- 
ing is somewhat steeper than those observed for molecular 
clouds. However this distribution of gas is an extreme case, 
as in general the gas would be comprised of multiple clumps 
of various sizes. We now consider many collisions and cal- 
culate the average velocity dispersion at each size scale. In 
all the collisions, we still assume 2 clumps collide, and fix 
the radius of the first clump to 1. However both the impact 
parameter b and the second clump radius, r2, can be varied. 
We show 3 different possibilities for the collisions of clumps 
in Fig. 17. Firstly we vary h randomly, assuming a proba- 
bility distribution function of /(6) = 1, < 6 < 1, with r2 
fixed at 1. Secondly, we take b biased toward larger values, 
so/(6)=6, 0<6<1, and r2 is again fixed at 1. Finally, we 
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Figure 17. The mass- weighted velocity dispersion predicted for 
the coUision of multiple spherical clumps. The solid line shows 
the velocity dispersion for 2 clumps of radii 1 and impact param- 
eter 0.5. The dot-dash and dotted lines give the average velocity 
dispersion for multiple collisions of clumps, where all clumps are 
assumed to have radii of 1, but the impact parameter varies. The 
impact parameter is either random (dot-dash) or biased towards 
large h (dotted). The velocity dispersion for multiple collisions 
with both random h and r2 is given by the long dashed line. The 
short dash line shows a oc r^-^ . 

allow both b and r2 to vary randomly, so f{b) = 1, < 6 < 1 
and /(r2) = 1, < r2 < 1. For comparison, the velocity size- 
scale relation from 2 clumps colliding with impact parameter 
6 = 0.5 (and both of radii 1) is also included. 

The slope is shallower for multiple collisions of random 
6, giving cr oc r°"^. The velocity dispersion was also calcu- 
lated with h biased towards large values, reflecting the prob- 
ability of a collision relative to the impact parameter. This 
produces a slightly flatter slope, although it mainly increases 
the magnitude of the velocity dispersion. Finally, with both 
h and r2 allowed to vary randomly, the gradient decreases 
further to a oc r°'^. 

This is still an incomplete model compared to the simu- 
lations of clumpy gas previously described, with further pos- 
sible variations in the clump and shock geometry. Further- 
more, collisions of multiple clumps have not been included, 
i.e. where > 2 clumps collide simultaneously, or 2 clumps 
collide with each other before colliding with further mate- 
rial in the shock. The latter is evident in the simulations, 
where layers of clumps enter the shock and interact with 
each other. This is again likely to produce a shallower rela- 
tion compared with the calculations in this section. Overall, 
for the most structured distribution of 2 offset clumps collid- 
ing, the resulting velocity size-scale relation is steep, a oc r^. 
For the least structured distribution, i.e. uniform gas, the ve- 
locity size-scale relation of the shock is flat. In reality, the 
distribution of gas, and therefore the gradient of the velocity 
size-scale relation, is likely to lie within these two extremes. 



5 CONCLUSION 

We have presented simulations and analysis of shocks for 
initial distributions of uniform, clumpy and fractal gas. We 
find an increasing velocity size scale relation in all our re- 
sults similar to that observed, except for the uniform shocks. 
For example, for a 2.2 D fractal distribution representative 



of the ISM, the velocity size-scale relation is approximately 
cr oc r°'^, in good agreement with observations. The slope of 
the velocity size-scale relation tends to decrease with distri- 
butions corresponding to higher filling factors. The velocity 
size-scale relations determined for these distributions can be 
understood in terms of mass loading, as indicated by our an- 
alytical tests. Oblique shock tests show that the magnitude 
of the Une-of sight velocity dispersion depends on the angle 
from which the shock is viewed, providing an observational 
test for this model. 

These results imply that: 1) The observed multi-scale 
structure of the ISM may explain the velocity dispersion in 
molecular clouds. This is in contrast to the usual view that 
turbulence produces the structure of molecular clouds (e.g. 
Falga rone et al.l (|2005D); 2) In these models, the dynamics 
of the shocked gas corresponds to random velocities, rather 
than classical turbulence, despite the apparent velocity size 
scaling relation. 
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